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Single-site and cluster dynamical mean-field methods are used to estimate the response of a doped 
Mott insulator to a charged impurity. The effect of correlations on the Thomas-Fermi screening 
properties is determined. The charge density, the on-site and near-neighbor spin-spin correlations 
in the vicinity of the impurity are compared to those far from it. The theory is used to address the 
question of the effect of the density perturbation induced by the muon charge on the local response 
functions of a high-temperature superconductor. For reasonable values of the background dielectric 
constant and basic correlation strength, a muon is shown to lead to an observable perturbation of the 
local spin dynamics, raising questions about the interpretation of muon-spin-rotation experiments 
in metallic high-temperature superconductors. 
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I. INTRODUCTION: 

Mobile electrons act to screen a charged impurity. 
Screening may be understood in terms of two equations: 
the Poisson equation which relates the electric poten- 
tial V(r) to the combination of the impurity charge den- 
sity Pimp('") and the change dn{r) in free charge density, 
and a constitutive equation which relates dn to V. For 
weakly correlated metals it suffices to linearize the con- 
stitutive relation so the screening properties are deter- 
mined by the density-density correlation function %. A 
locality approximation is typically appropriate, so that 
dn{r) = x(r = 0,uj = 0)V{r) and it is also reasonable 
to write the Poisson equation in its continuum form. 
These approximations imply that screening in weakly 
correlated metals is described by the familiar Thomas- 
Fermi equations, which lead to an exponential decay of 
the charge density characterized by the Thomas-Fermi 
length \tf = l/v/47re2x(r 0,w)/e. 

Remarkably, while the ubiquity of impurity effects in 
correlated electron materials has prompted extensive the- 
oretical studies of the consequences of local disorder for 
magnetic and superconducting properties,'^ few results 
seem to be available for the problem of screening of a 
charge center in a material with strong electronic corre- 
lations. Several effects appear to be important. First, 
and most trivially, lattice effects are strong so a discrete 
version of the Poisson equation must be used. Second, 
most correlated materials of interest are oxides or organ- 
ics with high background polarizability. Third, the con- 
stitutive relation between potential and density is likely 
to be strongly affected by correlation phenomena, which, 
in particular, will act to reduce the charge response. 
Fourth, in correlated materials, properties are typically 
sensitive functions of density, so that linearization in the 
magnitude of the density change may not be appropriate, 
while charge accumulation or depletion near an impurity 
may change the physics locally, for example, nucleating 
or suppressing local order or fluctuations. 

The possibility of local changes in the physics is of par- 



ticular importance in the context of muon-spin-rotation 
spectroscopy. In this class of experiments, a positively 
charged muon with a known initial spin direction is in- 
jected into a solid. Coupling to magnetic order or fluc- 
tuations causes the spin of the muon to precess before 
it decays and the amount of precession (and hence some 
information about the spin fluctuations) can be inferred 
from the angular distribution of the muon decay prod- 
ucts. If the charge of the muon causes a significant per- 
turbation of the electronic properties near the muon site, 
then the muon does not necessarily measure the intrinsic 
magnetic dynamics of the material. 

The question of the perturbation imposed by an in- 
jected muon has recently arisen in the context of the 
possible observation of an "orbital current" phase in 
high-temperat ure superconductors. Following a predic- 
tion of Varmajf^ neutron-scattering experimentgBEl re- 
ported evidence of a time-reversal symmetry-breaking 
phase characterized by local magnetic fields which are 
non-vanishing but average to zero over a unit cell; how- 
ever, muon-spin-rotation experiment^^^Sl failed to detect 
the magnetic fields implied by the neutron experiments. 
One possible resolution of the discrepancy is that the neu- 
tron measurements detect properties of a minority phase. 
Another possible resolution, proposed by Shekhter et 
al.^ is that the muon, which carries unit charge, per- 
turbs the local physics strongly enough to destroy the 
local order detected by neutrons. Reference [TO] pre- 
sented a Thomas-Fermi calculation which used a con- 
tinuum version of the Poisson equation, a value e ^ 4 
of the dielectric constant rather smaller than the value 
e ~ 10 — 15 generally accepted for oxides and a com- 
pressibility which was assumed to be linearizable and un- 
renormalized by many-body effects (although some con- 
sequences of the correlations were mentioned) . The cal- 
culation of Shekhter et al. implies that the muon would 
constitute a strong perturbation, dramatically changing 
the doping and the magnetic dynamics. However, the 
discussion given above implies that the assumptions on 
which the calculation is based may be questioned. 



In this paper we re-examine the issue. We use a tight- 
binding model description which captures the physics as- 
sociated with the discreteness of the lattice, we exam- 
ine the dependence on background dielectric constant, 
and most importantly we use single-sit^^ and cluster^ 
dynamical mean-field-based methods to provide an es- 
timate of the correlation effects on screening and on 
near-impurity electronic properties. We determine when 
linearization is appropriate and, where needed, use the 
full nonlinear (but local) charge response. We compute 
locally defined quantities which give some insight into 
changes in spin correlation properties. While our specific 
results are obtained for a lattice, doping, and interaction 
strength appropriate to high- Tc cuprates, we expect that 
our methods are more broadly applicable and our quali- 
tative results are relevant to a wider range of systems. 

We find that if dielectric constants in the physically 
reasonable range are used, the presence of a unit charge 
induces density changes which are a non-negligible frac- 
tion of the doping; however, the resulting changes in lo- 
cal magnetic properties are found to be modest, although 
observable. Further theoretical attention, perhaps using 
"LDA-KDMFT" mcthodJ^l^ should be given to model- 
ing the effects of charge centers in general and muons in 
particular. 
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FIG. 1: (Color online) Sketch of two cases considered in paper. 
Left panel: charged impurity (muon) placed at center of cube 
of transition-metal (Cu) sites. Right panel: impurity placed 
in plane. 



position. Placing the impurity at these high symmetry 
points allows us to use existing codes but as will be seen 
the physics we find is generic. 

We treat the screening using the self-consistent Hartree 
approximation. The presence of the impurity potential 
changes the electronic density on site i (coordinate Ri) 
from the average value n to a new value ni — n + brii s,o 
the total electrostatic potential is 



II. MODEL 

We approximate the conduction band degrees of free- 
dom as a one-band Hubbard model 

ij (7 % i 

(1) 

with hopping tij , on-site repulsion U and a spatially vary- 
ing electrochemical potential fii = p. + Vi determined self- 
consistently (see below) from the impurity charge and 
any induced electronic charge. For the explicit calcu- 
lations presented in this paper we take a set of planes 
which are electronically decoupled (no interplane hop- 
ping) but coupled via the Coulomb interaction. Each 
plane is taken to be a two-dimensional square lattice with 
nearest-neighbor hopping. We consider interactions of 
the order of the critical value Uc2 needed to drive a metal 
insulator transition in a homogeneous bulk system with 
one electron per site. /I, the chemical potential far from 
the impurity site, is chosen to produce a carrier density 
n corresponding to a hole doping oiS = l — n = 0.1 cor- 
responding approximately to the doping level at which 
pseudogap and magnetic effects occur in the high- Tc 
cuprates. 

We suppose that the impurity is located at a position 
Rfj, and has a charge e. We consider in detail two cases, 
shown in Fig. [T] In one we take the impurity to be lo- 
cated in the center of a plaquette in a Cu02 plane (i.e. 
the (1/2,1/2,0) position). In the other we place the impu- 
rity symmetrically between planes at the (1/2, 1/2, 1/2) 



v^ = — ^ — ^ + — 2J ^ ^ (2) 
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with e the background dielectric constant. The appropri- 
ate value of e is not well established. Optical conductiv- 
ity measurements-^'* suggest that e(w — > 0, g = 0) « 4, 
a value used in Ref. [101 however what is needed is 
e{uj = 0,q) for a range of q of the order of the shortest 
distance from a lattice site to the impurity position. Rea- 
sonable values of this quantity have not been determined. 
Calculations of the "screened U" for the related oxide ma- 
terial SrVOa yield a high frequency unscreened U ^ 14eV 
and a low frequency screened quantity W ~ 2eV suggest- 
ing an electronic contribution to e of e « 7. Lattice relax- 
ation effects may increase the short scale e to a number 
of order 15 (see, e.g. Ref. '15]) but of course lattice relax- 
ations may induce other changes in the model. Resolving 
these uncertainties is beyond the scope of the paper; we 
have therefore performed our explicit calculations for the 
two values e = 4 and 15. 

The remaining issue is the computation of 6ni for a 
given distribution of Vf, this is discussed in the next sec- 
tion. Here we note that the scale of the screening effects 
is set by the dimensionless parameter 

dn 

1=—T~ ^) 
ea d/i 

with a the in-plane lattice parameter. We use the value 
a « S.SA appropriate for cuprates. The band theory 
estimate for the compressibility dn/dfj, of weakly corre- 
lated electrons in the cuprate band structure is dn/dfi w 
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FIG. 2: (Color online) Sketch of the self-consistency proce- 
dure used to calculate the charge density and hybridization 
functions in the vicinity of a charged impurity. Starting from 
an initial guess for the potential Vi and the hybridization func- 
tion Ai{u}) the dynamical mean-field procedure is used to ob- 
tain a new density and hybridization function; the new density 
is used in Eq. ^ to obtain new potentials, and the process 
is iterated to self-consistency. Changes in the hybridization 
function are found to be sufficiently small that the DMFT 
loop may be solved once for the bulk material and dn/dfi 
obtained from this solution may then be used to update the 
density. 



1.4/eV SO that 7nonint ~ 5/e. As we shall see, for cor- 
relation strengths of the order of those believed to be 
relevant for high-temperature superconductors, the ac- 
tual compressibility, and therefore the actual 7 are likely 
to be about an order of magnitude smaller. 



III. METHOD 

We require the solution of a correlated electron prob- 
lem in a spatially inhomogcneous, self-consistently deter- 
mined potential. There is no general and exact method 
for obtaining this information. We adopt here the 
single-sit^S] and cellulai^i^ dynamical mean-field theory 
(CDMFT) approximations. These methods capture im- 
portant aspects of the strong correlation problem and, 
in particular, produce a Mott transition. The single-site 
method is more computationally tractable; hovifever the 
cluster method includes intersite correlations and may 
provide a more reasonable picture of the spin dynamics. 

In the single-site DMFT method, the electron self en- 
ergy E (in general a function of two coordinates and a 
frequency) is taken to be site local but may be different 
from site to site: E(i,j;a;) — )■ T,i{ui). The self energy on 
each site i is determined from the solution of a quantum 
impurity model (also different on each site). The impu- 
rity model is specified by the interaction Uri^n^ and a hy- 
bridization function Ai(w). The impurity model Green's 
function on site i is thus 

G'rH = [lo + /i + V,- A,(c.) - E,(c.)]-i . (4) 

The hybridization function is fixed by a self- 
consistency condition linking the Green's function 
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FIG. 3: (Color online) Local-density change induced by lo- 
cally potential In all cases bulk density is n = 0.9. 
Solid line: density change on site i, Srii, induced by poten- 
tial Sfii applied on same site, computed for non-interacting 
electrons using tight-binding band parameters appropriate to 
high-temperature superconductors. Dotted line: same com- 
putation, but using single-site dynamical mean-field theory 
with U = 13t and hybridization parameters taken from calcu- 
lation for U — 13t and 5 = 0.1. Dashed line: change in density 
on one site of a four-site plaquette induced by potential 5fi 
applied to all four sites of plaquette, computed using CDMFT 
dynamical mean-field theory with U = 9t with hybridization 
function corresponding to 17 = 9t and 5 — 0.1. 



G""'^(ti;) of the impurity model on site i to the ii compo- 
nent of the lattice Green's function, 



Gr(a;) = Gf (c.) 



I J i 



(5) 



There are now two issues of self-consistency: must 
be made self-consistent using Eqs. ^ and ([5| and the 
potentials Vi on all sites i must be made self consistent 
with the computed densities (which are obtained from 
the G™^(aj)) using Eq. To reach self-consistency one 
begins with an initial guess for the site densities. From 
this one computes the Vi via Eq. ([2| . Using these Vi and 
an initial guess for the hybridization function one solves 
the DMFT equations, obtaining converged solutions for 
G"""^ and A^. From these we recompute rii and hence V 
and continue the cycle until convergence is reached. 

Observe that the result of this procedure is that each 
site has a hybridization function determined by neigh- 
boring sites, which have different densities. Thus a given 
site "knows" that it is in a spatially inhomogeneous envi- 
ronment, and therefore has properties which are different 
from those of a hypothetical bulk system in which all sites 
have a density equal to the density of the designated site. 
At various points in the ensuing discussion we will com- 
pare properties of a given lattice site i with density to 
those that would be obtained in a bulk solid in which all 
sites had density n^. 

In practice the laborious procedure described above 
may be simplified. We have verified (for an example, see 
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Fig. 4]) that to within an accuracy of ~ 10% the change in 
the hybridization function is neghgible and the potential 
may be computed from the pure-system n(/Lt) curve. Thus 
we use the homogeneous bulk hybridization function to 
compute the variation in local n with local /i and use 
this to define the density/potential self-consistency with 
the result that it is necessary to solve the DMFT loop 
only once. The self-consistency procedure is sketched in 
Fig. (2). 

Figure [3] shows the change in on-site density caused by 
a change in local chemical potential computed for non- 
interacting electrons with bulk density corresponding to 
0.1 hole doping, and for two dynamical mean-field cases: 
single-site DMFT at U ^ lit (slightly larger than the 
single-site DMFT critical U which is \2t for this problem) 
and cluster DMFT at C/ = 9f which is rather larger than 
the Qt needed to open a gap in this approximation but is 
in the range believed to be reasonable for cuprates.^'' In 
the DMFT calculations the hybridization function was 
fixed at the form appropriate to a bulk material with 
density n = 0.9. As expected, the correlation effects sub- 
stantially reduce the local charge susceptibility: the ini- 
tial slope is decreased by a factor of about 5 relative to the 
noninteracting value and there is a substantial curvature. 
We also observe that even for 5n — 0.1, corresponding to 
a local density of one electron per site, the charge suscep- 
tibility remains non-vanishing (as expected because the 
local site is embedded in a metallic bath), whereas the 
corresponding bulk system with density n = \ per site 
would be in a gapped phase with vanishing compressibil- 
ity. 

The single-site method is reasonably computationally 
tractable, enabling the exploration of a wider parameter 
space and a relatively detailed computation of physical 
quantities. However, this approximation overestimates 
the critical interaction required to drive the Mott transi- 
tion, does not describe the "pseudogap" physics associ- 
ated with underdoped cuprates and more generally does 
not capture the physics associated with short-ranged in- 
tersite correlations. Cluster dynamical mean-field meth- 
ods capture more aspects of cuprate physics, including 
a lower critical value for the Mott transition and some 
aspects of intersite spin correlations. They also exhibit 
a pseudogap. However, the cluster methods are much 
more computationally expensive. Further, the widely 
used "dynamical cluster approximation"'^ is unsuited 
to impurity problems because it requires translation in- 
variance. We therefore adopt the CDMFT methocP® in 
which one tiles the lattice into real-space clusters; each 
cluster is regarded as a site of a new lattice of super- 
cells. The hopping terms connecting sites in the same 
supercell are part of the cluster Hamiltonian while the 
hopping terms connecting sites on different clusters de- 
fine the supercell band structure. The new lattice is 
treated via single-site dynamical mean-field theory (al- 
beit with a more complicated impurity), thus the self- 
consistency loop is the same as in the single-site case. 
We use four-site clusters. We solve the impurity model 



using the continuous-time quantum Monte Carlo method 
introduced in Ref. il8j for the four-site cluster we use 
the general ("matrix") representation of Ref. HH The 
method gives access both to the physical (lattice) elec- 
tron Green's functions and to correlation functions de- 
fined on the cluster model. While the cluster correlation 
functions are not identical to the corresponding lattice 
quantities, they are reasonable estimators of the physical 
correlators. 

One restriction should be noted: the impurity solver 
algorithm we us^^^ makes heavy use of symmetries and 
therefore requires that the four sites in the cluster have 
the same potential. Thus for cluster calculations we are 
limited to the case in which the impurity potential is 
the same for all four sites in the cluster. The geome- 
try we use guarantees that this is the case for the four 
sites closest to the impurity on each plane. However, for 
the farther plaquettes, a problem arises, because one side 
of a plaquette is necessarily closer to the impurity than 
the other, so the local symmetry is broken. We treat 
this situation by solving the Poisson equation and then 
on each cluster replacing the potential by the average of 
the calculated potential over the cluster sites. The long 
range of the Coulomb interaction and the relatively small 
changes induced on farther neighbor clusters make this a 
reasonable approximation. 



IV. RESULTS: DENSITY DISTRIBUTION 

An electrical conductor responds to a charged impurity 
by producing an electron density modulation ( "screening 
cloud") which screens the impurity charge. The panels 
of Fig. |4] present the spatial distribution of the screening 
cloud induced by an impurity of charge -t-1 in a hole- 
doped superconductor. Shown is the charge density per 
lattice site along a line passing near to the impurity site 
for two choices of background dielectric constant, e — A 
and e = 15, for two choices of impurity position (between 
planes or in the center of a plaquette in one plane) and for 
the two approximations we have used. We see that in all 
cases the density change is only appreciable on the sites 
adjacent to the impurity. For e = 4 the density change 
on the sites nearest to the impurity is large enough to 
move the local density very close to the half filled value. 
For e = 15 the density change is about a factor of two 
smaller than for e = 4. The density profiles calculated for 
single-site and four-site DMFT are very similar, because 
the density profile is controlled by the local compressibil- 
ity which is similar for the two cases we have considered. 
The density profiles calculated for the two impurity lo- 
cations and for the farther plane are also similar because 
the 1 /r variation in the unscreened Coulomb interaction 
is relatively slow. These calculations are performed using 
the simplified self-consistency loop described above; also 
shown are results obtained using the full DMFT proce- 
dure for the e = 15, single-site DMFT case. The density 
changes induced by the impurity potential are seen to 
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FIG. 4: (Color online) (a-d) Variation in conduction-band 
density per lattice site along line {x, 0, 0) induced by impuri- 
ties positioned at [(a) and (c)] (1/2,1/2,1/2) and [(b) and 
(d)] (1/2,1/2,0) calculated as described in the text using 
single-site DMFT with U = 13t and and four-site DMFT 
with U = 9t. [(e) and (f)] Variation in conduction-band den- 
sity per lattice site along line {x, 0, 1) induced by impurity 
positioned at (1/2,1/2,0) for (e) e = 4 and (f) e = 15 cal- 
culated as described in the text using single-site DMFT with 
U = 13t and four-site DMFT with U = 9t. 



be generically of the order of 0.05 electrons per site or 
less, which is less than but of the order of the doping for 
underdoped high-temperature superconductors. 



V. RESULTS: SPIN CORRELATIONS 

In this section we study how the screening cloud affects 
the local spin dynamics. This is not straightforward be- 
cause the spin dynamics are expected to be strongly dop- 
ing dependent in a homogeneous bulk system while here 
we must treat a spatially inhomogeneous system. We 
study impurity-model correlation functions, which can be 
directly measured in our simulations. These are not iden- 
tical to the spin correlations of the actual lattice problem 
but are expected to have similar magnitude and similar 
doping and temperature dependence to those of the full 
lattice problem. Further, the simulation gives results for 
Matsubara frequencies w„ = 2'KnT. The n = term is 
in essence the classical (thermal) part of the spin-spin 
correlation function while the n 7^ terms give some 
information on the quantum fluctuations in the system. 
We present results for the spin correlations on the site 
nearest to the charge center and for the second neigh- 
bor, and compare the results to those found far from the 
charge center and also to those computed for a hypo- 
thetical bulk system with average density equal to that 
on the site nearest to the charge center. 
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FIG. 5: (Color online) Impurity-model spin-spin correlation 
along line (x, 0, 0) computed from single-site dynamical mean- 
field theory as function of Matsubara frequency at inverse 
temperature 13 — 20 and dielectric constants indicated for 
charged impurity at position (1/2, 1/2,0). Upper panels are 
onsite correlators for sites nearest to muon {x = 1), lower 
panels are onsite correlators for second neighbors to muon 
(a; = 2), correlator for site far from muon is also included. 



Figure [5] shows the impurity-model spin-correlation 
functions computed using single-site dynamical mean- 
field theory. We see that the "classical" (zero Matsubara 
frequency) spin correlations of sites near the muon are 
enhanced relative to the value far from the muon site 
but are not as large as those of a hypothetical system 
with average density equal to the density on the near- 
muon sites. This shows that the spin correlations on a 
given site are controlled not only by the density on the 
site but also by the properties of the neighboring sites. 
For the case e = 4 we see that the changes are substantial 
(increasing the value at the lowest Matsubara frequency 
by a factor of about 1.5, which in turn is about half of 
the increase that would occur in a sample whose aver- 
age density was set equal to the density on the impurity 
site). On the other hand for e = 15 the changes, although 
visible, are much smaller. 

The single-site dynamical mean-field theory is known 
to provide a poor approximation to the spin correlations 
of a doped Mott insulator, at least in two spatial dimen- 
sions, because it neglects antiferromagnetic correlations. 
We have therefore also considered the spin correlations 
in the four-site CDMFT calculations. Here the differ- 
ence between on-site and first-neighbor impurity-model 
spin correlations reveals the importance of antiferromag- 
netic fluctuations. Figure [6] shows results for the same 
parameters as in Figure [5] We see that the first-neighbor 
correlations increase by about 40% for e = 4 and about 
25% for e ~ 15. We have also examined other correla- 
tion functions, in particular the equal-time singlet-pair 
correlations which are the dominant fluctuations on the 
four-site plaquette, finding that these are enhanced by 
similar amounts. 
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FIG. 6: (Color online) On-site (upper panels) and first- 
neighbor (lower panels) impurity model spin-spin correlation 
computed from four-site CDMFT cluster dynamical mean- 
field theory as function of Matsubara frequency at inverse 
temperature /3 = 20 and dielectric constants indicated, for 
charged impurity at position (1/2, 1/2, 0). Solid line with solid 
squares (blue online) shows correlator computed for sites far 
from impurity site. Dotted line with solid circles (red online) 
shows correlator computed on designated site. Dashed line 
with open squares (blue online) shows correlated computed 
for bulk system with mean density equal to density on impu- 
rity site. 



VI. CONCLUSION 

In this work, we have shown how "strong correlation" 
effects alter the response of a material to a local charge in- 
homogeneity. We introduced a general method, based on 
dynamical mean-field theory, for calculating these effects 
and applied it to the question of the changes produced by 
the presence of a muon in a high-temperature supercon- 
ductor. We found, in qualitative agreement with previous 
work,^° that the muon is not a "soft probe:" although 
the main correlation effect is a suppression of the charge 
susceptibility by a factor of 3-4 relative to band theory, 
the charge field associated with the muon may produce a 
significant change in the charge density on nearby sites, of 
order 0.05 electrons per site. In a less strongly correlated 
material, the change in charge density on the near-muon 
sites would be larger. A crucial issue in determining the 
scale of the effects was found to be the value of the static, 
short wavelength, dielectric constant. The effect on the 
local spin correlations is smaller than the effect on the 
charge density but is not negligible. 

Our calculation involves several approximations. The 
most important is the value of the dielectric constant, 
which is not known a priori. Varying the dielectric con- 
stant over the range of values which have been proposed 
for cuprates or other oxides leads to factor of 2 changes 
in our results. A closely related issue concerns the local 
lattice distortions which would normally be induced near 
a charged impurity (see Ref. [15] for an example in a dif- 
ferent context) . These have the potential to change local 
hopping amplitudes and perhaps correlation strengths, 
although correlation strengths, being an atomic property. 



will be less strongly affected. In view of the importance 
of muon-spin rotation as a probe of correlated materials, 
these issues deserve further investigation, perhaps via a 
band-theory calculations. A second approximation is the 
use of the dynamical mean-field method. The consistency 
of our single-site and four-site results for the charge com- 
pressibility suggest that our basic findings for the den- 
sity correlations are a reasonable estimate of the correct 
behavior. However, the calculated spin correlations are 
probably subject to larger uncertainties, which are at this 
point not easily quantified. We know that the charge per- 
turbation is important only on sites immediately adjacent 
to the charge center. The spin correlations on these sites, 
which are the ones which would be probed by a muon, 
are affected both by the on-site property (the change in 
local density) and by the properties of the nearby sites, 
and the nearby sites in turn both affect and are affected 
by the near-impurity sites. If the intersite spin coupling 
is strong it is possible that the magnetic properties are 
controlled by the sites farther away from the muon. It 
is very likely that the dynamical mean-field methods we 
use underestimate these spin-correlation effects. Their 
investigation is an important open problem. 

While our specific numerical results were obtained for 
model parameters appropriate to high-temperature su- 
perconductors, they have implications for the more gen- 
eral issue of the response of a correlated electron mate- 
rial to a charged impurity. To illustrate this point we 
consider a generic correlated material, which we assume 
to be a more or less cubic lattice. For simplicity place 
the charged impurity at the center of a cube of sites. 
The charged impurity will induce a screening cloud con- 
taining one electron. The length scale over which this 
charge is distributed is set by the density-density cor- 
relation function of the charged material and the back- 
ground dielectric function. If values typical of a weakly 
correlated material are used and the dielectric function is 
of the order of 10 or less, a simple extension of the esti- 
mates we have presented indicates that almost the entire 
screening charge sits as close to the impurity as it can get. 
The density change on the near-impurity sites would then 
(for the simple cubic situation we have considered) be ap- 
proximately 1/6 electron per site, a change large enough 
to affect the local physics. On the other hand, if strong 
correlation effects are important (as in the case of high- 
Tc materials where they reduce the charge response by 
a factor of 5 or more) , the total charge would be spread 
over a wider range and the concentration on the near- 
impurity sites would be substantially smaller. However, 
the relative effect on the local physics would still not 
be small, as the greatest suppression of charge response 
occurs for a lightly doped Mott insulator, where the im- 
portant scale is the doping, which would itself be small. 
Thus even in this case we would expect that a charged 
impurity would change the local physics noticeably. A 
quantitative test of our theory would involve measure- 
ments of the near-impurity charge-density profile and a 
comparative measurement of spin dynamics near to and 



7 



far from the impurity site. 

Our results have impUcations for muon-spin-resonance 
experiments on transition-metal oxides. Muons are an 
important probe of the spin dynamical of condensed mat- 
ter physics, but a muon has a charge -1-1, and the re- 
sults presented here indicate that in transition metal ox- 
ides a muon is unlikely to be a "soft" probe; rather, it 
significantly perturbs the medium in which it is embed- 
ded. Feyerherm et al^^ reached a similar conclusion in 
a study of PrNis, a rare earth system with more compli- 
cated physics, showing that muons significantly perturb 
the crystal-field structure on the near muon Pr sites. In 
the case studied here we showed that the perturbation 
due to the muon affects local properties such as the near- 
muon spin dynamics and presumably (although we have 
not investigated this) the size of the ordered moment. 
It is important to note that a dilute concentration of 



muons should have only negligible effects on "global" or 
long-range properties such as magnetic phase boundaries 
or superconducting penetration depth. 
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